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A NEW METHOD EOR CONSTRAINING MOLECULAR CLOUD THICKNESS: A STUDY OE TAURUS, 

PERSEUS AND OPHIUCHUS 

Lei Qian \ Di Li ^ ^ '‘j Stella Offner ^ and Zhichen Pan ^ 


ABSTRACT 

The core velocity dispersion (CVD) is a potentially useful tool for studying the turbulent velocity 
field of molecular clouds. CVD is based on centroid velocities of dense gas clumps, thus is less prone to 
density fluctuation and reflects more directly the cloud velocity field. Prior work demonstrated that the 
Taurus molecular cloud CVD resembles the well-known Larson’s linewidth-size relation of molecular 
clouds. In this work, we studied the dependence of the CVD on the line-of-sight thickness of molecular 
clouds, a quantity which cannot be measured by direct means. We produced a simple statistical model 
of cores within clouds and analyzed the CVD of a variety of hydrodynamical simulations. We show 
that the relation between the CVD and the 2D projected separation of cores (L2 d) is sensitive to the 
cloud thickness. When the cloud is thin, the index of CVD-L2ri relation (7 in the relation CVD~ 
reflects the underlying energy spectrum {E{k) ^ k~^) in that ^^{13 — l)/2. The CVD-L2L1 relation 
becomes flatter (7 —0) for thicker clouds. We used this result to constrain the thicknesses of Taurus, 
Perseus, and Ophiuchus. We conclude that Taurus has a ratio of cloud depth to cloud length smaller 
than about I/I0-I/8, i.e. it is a sheet. A simple geometric model fit to the linewidth-size relation 
indicates that the Taurus cloud has a ^ 0.7 pc line-of-sight dimension. In contrast, Perseus and 
Ophiuchus are thicker and have ratios of cloud depth to cloud length larger than about 1/10-1/8. 

Subject headings: ISM: clouds — ISM: molecules 


1. INTRODUCTION 

Stars form in molecular clouds. There is much evi¬ 
dence that molecular clouds are turbulent, e.g. the ob¬ 
served molecular linewidth is much larger than the ther¬ 
mal linewidth. Because the interstellar medium (ISM) 
is turbulent, its velocity spectrum is self-similar, i.e., a 
power-l aw, as evidenced in numerous studies (see the re¬ 
view bv iHennebelle fc Ealgaroii3l2012L and the references 
therein). 

Turbulence in molecular clouds has be e n stu died in 
a variety of scales and contexts. iLarsonl (jI98Hl found 
for the first time a power law relation between the ve¬ 
locity dispersion a and the (projected, 2D) size L of 
molecular clouds oc This power law form re¬ 

lation between the velocity dispersion and the size of 
molecular clouds is then called Larson’s relation. Later 
observation s established the now widely accepted in¬ 
dex of 0.5 (iSolomon et all Il987t lEalerarone et al.l 119921 : 
iHever fc BrnntJ 120041 . The former index is close to the 
power law index 1/3 of incompressible turbulence, while 
the latter m ay indicate density fluctuations wit hin molec¬ 
ular clouds (|Vazauez-Semadeni fc Gazollll995[l. which is 
confirmed in observational studies ( Brunt &: Heved[2002l : 
iBrnntI 1200311 and numerical simulations with increas¬ 
ing sophisticated physic s and resolution (lOstriker et al.l 
120011 lOffner et al.ll2008t iKritsuk et al.ll2013l l. However, 
recent studies on this topic bring more complications. 
iZhang et alJ ()2014[) find in North American and Pelican 
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molecular clouds a power law index of 0.43 for the Lar¬ 
son’s relation, which de viates from the trad itional value 
of 0.5. In another study, iHever et al.l (1200911 showed that 
the normalization of the velocity dispersion-size relation 
molecular clouds, uq = , scales with the surface 

density of the clouds S as uq ~ indicating that 

molecular clouds are in self-gravitational equilibrium. 

Traditionally, str ucture function methods (e.g. 
iBrunt fc Hevedl2002ll and princ i pal compone nt analysis 
(PCA. IHever fc Schloerbl Il997t iBruntI 1200311 are used 
for charact e rizing the turbulent velocity field. Recently, 
iQian et HI (j2012[ l developed the core velocity dispersion 
(CVD) technique to study turbulence in the Taurus 
molecular cloud. This method uses the centroid ve¬ 
locity (average velocity of a core weighted with mass) 
of identified molecular cores to sample the velocity 
field of spatially and spectrally coherent structures in 
relatively dense gas. In this paper, we compare the 
characteristics of the CVD with the gas turbulence 
and examine whether CVD could trace the general 
motion of the underlying gas. IQian et ^ (|2012h found 
that the dispersion of the core velocities increases with 
distance, i.e. CVDoc between 1 pc and 10 pc, which 
resembles Larson’s relation (hereafter, we will attach 
the subscript 2D to L as L^d and 3D to I as I^d to 
emphasize the dimension of the spatial scales). Larson’s 
relation can be explained b y the velocity spectr um of a 
3D turbulent velocity field (|Kritsuk et al.l[2M^ . There 
is a fundamental difference between the CVD-L2D 
relation and the turbulent velocity spectrum, in that the 
turbulent velocity spectrum is three-dimensional, while 
the spatial scale in the CVD-L2 _d relation is a projected 
(2D) scale. 

In this work, we expanded upon the previous results 
and study the influence of the line-of-sight scale (thick- 
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ness) on the CVD-L 2 _d relation. We described our ob¬ 
servational data in presented our methods in Sj3l and 
discussed results in 2] In ^ we summarized our con¬ 
clusions and addressed the potential application of the 
CVD in future studies. 

2. OBSERVATIONAL DATA 

The (J=l-0, 110.2014 GHz) data of Taurus 

in this work were obtained with the 13.7 m FCRAO 
telescope between 2003 and 2005. The map is cen¬ 
tered at RA(2000.0) = 04^32'"44.6®, Dec(2000.0) = 
24°25'13.08”, with an area of ^ 98 deg^, and a noise 
level of 0.1 K (jNaravanan et all 1200811 . The velocity 
resolution is 0.266 km/s. The data of Perseus 

and Ophiuchus ca me from the COMPLETE database 
(|Ridge et al.ll2006lFI . which were also observed with the 
FCRAO telescope. The velocity resolution is 0.066 km/s. 
The noise levels of Perseus and Ophiuchus are 0.15 K and 
0.2 K, respe ctively. 

Following lOian et ^ (|2012f) . we used the GAUSS- 
CLUMPS method in the Starlink software^ to identify 
cores from the data cube. The lower thresholds of the 
peak inte nsity of cores we re set to be 7 times the noise 
level Isee lOian et al.l[2012[ l. which were about 0.7 K, 1.0 
K, and 1.5 K for Taurus, Perseus, and Ophiuchus, re¬ 
spectively. In Taurus, Perseus, and Ophiuchus, 588, 693, 
and 270 cores were identified, respectively. The typical 
size of the cores is around 0.1 pc. 

3. METHODS 

In order to investigate the influence of the cloud thick¬ 
ness on the CVD-L 2 _d relation, we used two approaches. 
In the first approach we statistically generate a sample 
of cores in a velocity field with second order structure 
function of index 0.5, which corresponds to an energy 
spectrum of index /? = 2. In the second approach we 
perform hydrodynamic simulations of molecular clouds. 

3.1. Calculation of the CVD 

CVD stands for Core Velocity Di spersion, which is ob¬ 
tained as follows (|Qian et al.ll2012ll . First the difference 
of centroid velocity (Sv) and projected spatial distance 
{L 2 d) of each pair of cores are calculated. The distribu¬ 
tion of Sv in the Sv — L 2 D plane is in grey scale in the top 
panel of Figure [T] for the Taurus molecular cloud. Then 
we calculate the CVD, which is the mean square root of 
the values of 5v, i.e. CVD =< 5v^ at different L 2 D 
(see the green diamonds in the top panel of Figur e ID)- 
Figured] is similar to Figure 18 in I Qian et al.l (I2012f) . but 
now we focus on the power-law-form CVD-L 2 D relation 
between 1 pc and 10 pc, which reflects the pattern of 
turbulent motion in Taurus. 

In order to compare the CVD to a turbulence velocity 
spectrum, we fit a power law or a broken power law to 
the CVD-L 2 _d relation in this work. The break point of 
the broken power law is determined from the log (CVD)- 
log(L 2 _D) plot by eye. The fitting error i s estimated with 
a bootstrap method ([Press et al.! Il992li . We construct 
100 samples of core pairs based on the real sample of 
core pairs, in which every pair in each sample is randomly 

® http://www.cfa.harvard.edu/COMPLETE/ 

® http://starlink.jach.hawaii.edu/starlink/ 


selected from the data. Then we fit parameters to these 
100 samples, and the error is estimated by calculating the 
root of mean square of the fitted parameters. Since there 
is degeneracy between the normalization and the index 
of a power law function, for power law fitting through 
out this paper, the power law index is first obtained by 
fitting a linear function to the log-log data and then fixed 
in the subsequent fitting in linear space, meant to obtain 
the offset constant. 
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Figure 1. Top: Core velocity difference between each core pair, 
5v^ and the core velocity dispersion (CVD= vs. the pro¬ 

jected separation, 1/21), of the cores for the Taurus molecular cloud. 
The background is an image of the number of data points in the 
5v — L 2 D plane, with the grey scale (density bar at the top of the 
figure) showing the density of points. The green diamonds repre¬ 
sent the CVD in each separation bin. The CVD data points with 
apparent separation 0 < L 2 D ^ 10 pc can be fitted with a power 
law of CVD (km/s)= (0.3 zb 0 . 04 )L 2 _d( pc)0-5±0-05 + ( 0.4 =b 0.05). 
The horizontal line shows the mean CVD value, 1.1 km/s, of data 
points with L 2 D > 10 pc. Bottom: The CVD as a function of pro¬ 
jected distance. The range of the power law fit is determined in this 
log-log scale plot. The fitting error is estimated with a bootstrap 
method. 


3.2. Statistical Core Distribution Model 

In order to study how well the core velocities cor¬ 
relate with the underlying velocity field and the cloud 
geometry, we formulate a simple statistical model, the 
fractional Brownian motion model, of cores within a 
turbulent molecular cloud. The core sample is gener¬ 
ated with procedures similar to those used in gener¬ 
ating a turbulent velocity fi eld ([Brunt &: Heved 120021: 
iMivillc-Deschenes et al.|[2003[) . Theoretically, the root of 
mean square turbulent velocity at scale I (scale length in 
3D), vi , can be determined by the velocity fiuctuations 










































Constraining Molecular Cloud Thickness 


3 


over scales smaller than 1: 
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E{k)dk oc 


k ^dk oc 


12-k/I 


( 1 ) 


where k is the wave number. If we assume < vf P, 

the relation between /3 and 7 will be 7 = — l)/2 in 

3D. By definition, 


J v{k)‘^k'^dk , 


( 2 ) 


where v{k) oc k~^. 

We first generate a 3D velocity field and put roughly 
4000 cores on a randomly selected set of grid points. 
A simple way to do this follows the procedure below. 
First, we generate a Gaussian random field with dimen¬ 
sion 128x128x128 on a 3D grid. Next, we perform a 
Fourier transform. Then, we process this field in the 
frequency domain to satisfy the desired power law distri¬ 
bution (oc k~^). Finally, we perform an inverse Fourier 
transformation and normalize the generated field to ful¬ 
fill the desired variance, which is related to the velocity 
dispersion. 


3.3. Hydrodynamic Simulations and Synthetic 
Observations 

We use hydrodynamic simulations to explore the ef¬ 
fects of density/velocity correlations. We use the 
ORION adaptive mesh refinement (AMR) code to perform 
simulations of turbu l ent star-form ing molecular clouds 
(|Truelove et ahlfTOOSt lKleinlll999D . The simulations use 
periodic boundary conditions and are meant to represent 
a piece of a typical low-mass star forming cloud similar 
to Taurus. The numerical methods we use are simila r 
to previous calculations fe.g. lOffner et al.ll2009bl . [20131 1. 
so we briefly describe the simulations here and refer the 
reader to these papers for additional details. 

We compare the observations with snapshots from two 
simulations, each with different turbulent properties (see 
Table [Ij. Both simulations have a domain size of 10 pc, 
contain ^ 15,OOOM0, and have an initial gas temper¬ 
ature of 10 K. Simulation RT has a turbulent velocity 
dispersion of 1.52 km s“^, such that it satisfies the ob¬ 
served linewi dth-size relation at 10 pc , cr 2 D = 0.72i?°-^ 
km s“^ (e.g., iMcKee fc Ostrikeil 1200711 . RTl and RT2 
correspond to two different output times for the RT sim¬ 
ulation. They have ages of 0.65 Myr and 1.03 Myr, re¬ 
spectively. The turbulence is initialized by adding ran¬ 
dom velocity perturbations with an input power spec¬ 
trum P{k) oc for input wave numbers in the range 
fc = 1 2. Simulation HD has a velocity dispersion 

of 1.06 km s“^ and the turbulent perturbations have 
P{k) oc k~^ for fc = 1 10. The latter turbulent per¬ 

turbations are generated in the same way as those used 
in section o To achieve a turbulent steady-state, in 
which the power-spectrum and density distribution func¬ 
tion are constant in time, we inject the perturbations for 
two domain crossing times without self-gravity. 

After turning on gravity, we continue injecting energy 
so that the velocity dispersion remains constant. With¬ 
out continuous energy injection, the turbulence would 
decay as a result of dissipation in shocks ov er a crossing 
time le.g. iStone et al.l[l998t lMacLowlll999D . Once col¬ 
lapse begins, refinement is added to ensure that the Jeans 


criterion is satisfied f or a Jeans number of Nj = 0.125 
l|Truelove et al.l Il997f) . We insert sink particles when 
the density exceed s the Jeans resolution on the max¬ 
imum AMR level ([Krumholz et al.l I2004fl . HD adopts 
an isothermal equation of state, while RT solves the 
radiative transfer equation in the flux limited diffusion 
approximation; it a lso includes radiative f e edback from 
the fo rming stars (|Krumholz et al.l l2007t lOffner et al.l 
l2009bf) . This heating, which is generally limited to re¬ 
gions near sink particles, has minimal impact on the tur¬ 
bulent properties. 

RT and HD have 512^ and 256^ base grids with 2 and 
4 AMR levels of refinement, respectively. Here, we only 
use the data on the base grid, which is comparable to the 
observational pixel resolution, to produce the synthetic 
observations (see below). Figure [5] shows the velocity 
energy spectra for the three snapshots. For low wave 
numbers (fc = 1 ~ 4) the energy spectrum is flatter than 
fc“^. It steepens above fc ^ 20 where dissipation occurs. 
The dissipation scale is set by the numerical grid size and 
typically occurs on a scale of a few grid cells. 



Figure 2. Simulated velocity energy spectra compensated by P 
for the hydrodynamic simulations. The grey dotted line indicates 
the inertial range. 

In order to compare the simulations directly with the 
Taurus observations we post-process the outputs with the 
radiative transfer code RADMC- 3e1!|. We use the Large 
Velocity Gradient (LVG) approach, which computes the 
molecule level popu lations by solving for radiative sta¬ 
tistical equilibrium (jShettv et al.ll201lh . The input den¬ 
sities, temperatures, and velocities are supplied by the 
hydrodynamic simulations. 

We model unresolved turbulent broadening by adopt¬ 
ing a constant micro-turbulence parameter of 0.1 km s“^. 
The H 2 number density is defined as tihs = p/(2-8mp), 
where p is the gas density, the 2.8 factor accounts for the 
molecular weight of Helium. We assume a constant ^^GO 
abundance of [^^GO/H 2 ]= 2 .5 x 10“® (|Frerking et al.l 
1 19821 : [hanger k. Penziaslll99^ . and we adopt t he mo lecu- 
lar collisional coefficients from iSchoier et al.l (l2005fl . We 
compute the line emission for the velocities within ±5 km 

^ http://www.ita.uni-heidelberg.de/ dullemond/software/radmc- 
3d/ 
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s“^ of the line rest frequency and use a channel resolution 
of An = 0.039 km s“^ . 

To account for differences in resolution and sensitivity, 
we convolve the synthetic line data with a 45” Gaussian 
beam, assuming that the simulated cloud lies at a dis¬ 
tance of 140 pc. We re-bin the emission such that the 
pixel size is 20” and the velocity channel width is 0.26 
km s“^. Finally, we add Gaussian random noise with a 
standard deviation of 0.3 K. 

In order to assess the effects of cloud thickness on the 
observations, we perform the line transfer on input vol¬ 
umes of different extents along the line of sight. We com¬ 
pare the cases for line-of-sight thicknesses L, ^L, jL, ^L, 
and where L = 10 pc. Table [5] summarizes all the 
RADMC-3D runs. 

4. RESULTS 

4.1. The CVD of Taurus, Perseus, and Ophiuchus 

The observed Taurus CVD for cores identified in 
^^CO(l-O) is shown in Figure[T] In the scale range L 2 D < 
10 pc, the CVD-L 2 _d relation is well-fit by a power law 
given by CVD= (0.3±0.04)L2 d(pc)°' 5±° °'5-k (0.4±0.05) 
km s“^. The characteristic scale of 10 pc can be clearly 
seen from the log scale plot in the bottom panel of Fig¬ 
ure [TJ 

The solid line in Figure |3] shows the histogram of 
the core centroid velocities. The distribution is double- 
peaked, which could indicate two core populations. This 
might occur if the Taurus has distinct sub-regions with 
different line-of-sight velocities. One group has centroid 
velocities larger then 6 km/s and the other has centroid 
velocities smaller than 6 km/s. The former group of cores 
does not have a clear velocity gradient, but the latter 
group does (Vw ~ 1 km/s/deg). A plot of the core cen¬ 
troid velocity versus right ascension as shown in Figured 
also indicates that there are two core populations. 

We performed the GVD analysis of each group of cores 
and found that the CVD-L 2 L 1 relation for cores with cen¬ 
troid velocities larger than 6 km/s can be fitted by a bro¬ 
ken power law with power law indices 0.1±0.03 {L 2 D < 5 
pc) and 0.3 ±0.02 (5 pc< L 2 D < 10 pc). Cores with cen¬ 
troid velocities smaller than 6 km/s can be fitted with 
an index of 0.5±0.3 (Figure[5|). We fit and then subtract 
the velocity gradient 1 km/s/deg in the latter case and 
find that the resultant power law index is also 0.5. 

We have performed a similar analysis for the Perseus 
and Ophiuchus molecular clouds. In Perseus, the cores 
do not divide into groups (dash-dotted line in Figure [3]), 
although there is a velocity gradient of 4 km/s/deg 
among the cores (see Figure IHl). We perform the fit 
again after subtracting the gradient and verify that the 
gradient does not affect our results. The CVD between 
0 pc < L 2 D < 7 pc is nearly constant, about 1.5 km/s, 
while the upturn between 7 pc < L 2 D ± 13 pc can 
be fitted with CVD(km/s)= (0.3±0.001)L2 d(pc)° ®=^° °i 
(F igure [7]). The break point (7 pc) is determined from 
the lower panel by eye. In Ophiuchus, the core veloc¬ 
ity distribution is also centrally peaked (dotted line in 
Figure[3]). There is no core velocity gradient (FigureE]), 
and we find that the CVD between 0 pc < L 2 D ± 3.5 
pc is nearly constant, about 0.8 km/s. The upturn in 
the range of 3.5 < L 2 D < 5 pc can be fitted with 
CVD(km/s)= (0.5±0.05)L2i?(pc)°-^=^°°’’-(0.03±0.001) 



0 2 4 6 8 10 12 14 
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Figure 3. Solid ZmeiHistogram of the central velocities of cores 
in the Taurus molecular cloud. There are clearly two groups of 
cores with different central velocities, divided by the solid line at 
6 km/s. Dash-dotted line: Histogram of the central velocities of 
cores in Perseus. Dotted line: Histogram of the central velocities 
of cores in Ophiuchus. 
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Figure 4. The centroid velocity of cores vs. RA. relation in the 
Taurus molecular cloud. There are two apparent core groups with 
the dividing line roughly lying at 6 km/s (solid line). One group 
of cores has a velocity gradient. 

(Figure[Hl). The break point (3.5 pc) is determined from 
the lower panel by eye. 

The breaks in Perseus (7 pc) and in Ophiuchus (3.5 
pc) may indicate the thicknesses of the clouds, which 
are discussed in section 14.21 

4.2. Influence of geometry 

Previous studies of the Taurus molecular cloud found 
that the CVD follows a power law in the projected length 
scale (2D distanc e, L 2 d) with a power law index of 0.5 
(jOian et al.l 1201^ . It resembles the Larson’s relation. 
The main difference between these two relations is that 
the Larson’s relation is based on the integrated proper¬ 
ties of different clouds, while the centroid velocity of the 
cores used in the CVD calculation probes the systemic 
velocities of spatially and spectrally ’coherent’ clumps of 
gas. 

In principle, if a cloud has an infinite thickness (but 
finite optical depth), the CVD-L 2 D will not hold. The 
only reason that there is a power law relation between 
the CVD and the projected length scale L 2 D is that 
the cores are distributed in a limited range along the 
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Table 1 

Simulation Properties 


Model‘s 

M(Mq) 

M 

L (pc) 

t/(Myr) 

ATS 

Axniinfpc) 

RTl 

1.475 x 10"* 

14.0 

10.0 

0.65 

5123 

0.005 

RT2 

1.475 x 10^ 

14.0 

10.0 

1.03 

5123 

0.005 

HDl 

1.475 x 10^ 

10.0 

10.0 

1.06 

2563 

0.0025 


^ Run name, total gas mass, domain length, analysis output time, base- 
grid size, and minimum cell size. RTl and RT2 are two different output 
times for the same simulation. They have random perturbations in the 
range = 1 ~ 2 with P(fci_ 2 ) cc HDl has random perturbations 
fc = 1 10 with -P(Aii^io)) oc k~^. All runs are first evolved for two 

crossing times without self-gravity and have an initial temperature of 
10 K. RTl and RT2 also include radiative feedback from forming stars 
(Offner et al. 2009). 


Table 2 
RADMC Runs 


ModeD 


X 

Ny 

X 

N, 

View 

(Lxipc), Ly{pc), Lzfpc)) 

RTl_L10z_M14 

512 

X 

512 

X 

512 

z 

(10, 10, 10) 

RTl_L5z_M14 

512 

X 

512 

X 

256 

z 

(10, 10, 5) 

RTl_L2.5z_M14 

512 

X 

512 

X 

128 

z 

(10, 10, 2.5) 

RTl_L0.125z_M14 

512 

X 

512 

X 

64 

z 

(10, 10, 0.125) 

RTl_L0.0625z_M14 

512 

X 

512 

X 

32 

z 

(10, 10, 0.0625) 

RTl_L10y_M14 

512 

X 

512 

X 

512 

y 

(10, 10, 10) 

RTl_L5y_M14 

512 

X 

512 

X 

256 

y 

(10, 10, 5) 

RTl_L2.5y.M14 

512 

X 

512 

X 

128 

y 

(10, 10, 2.5) 

RT2_L10y.M14 

512 

X 

512 

X 

512 

y 

(10, 10, 10) 

HDl.L10z.M10 

256 

X 

256 

X 

256 

Z 

(10, 10, 10) 

HDl.L5z.M10 

256 

X 

256 

X 

128 

z 

(10, 10, 5) 


^ Model name, input RADMC grid size, line-of-sight view, and domain size. 
We adopt a constant micro-turbulence value of 0.1 km/s and a doppler 
catching parameter of 0.25 for all runs. 


line of sight direction (’thickness’), in which case 2D 
still bears resemblance to 3D. The difference between 
2D and 3D description of turbulence have been studied 
in som e works based on the line centroid velocity analy¬ 
sis (e .g. lMiville-Deschenes et al.l[2003HBrunt &: Mac Lowl 
I2004f) . In the case where the line of sight scale is much 
smaller than the transverse scale, the exponent 72 D 
of the relation between centroid velocity disper¬ 
sion (ay) and projected scale L 2 D is related to the 
exponent o f the second order structure f unction by 
72 D = 7313 (iMiville-Deschenes et al.l[200l) . Here, 
is the power law index of the <Jv-L 2 d/ 3 D relation (section 
IQ) . However, in general, the density inhomogeneity will 
produce 7213 « 7313 « 0.5. When the effect of the den¬ 
sity inhomogeneity is diminished in the late time of de¬ 
caying turbulence, the relation between the 2D power 
law index 720 and the S D one 73 D is 7213 ~ -\- 0.5 

(|Brunt &: Mac Lo^ l2004fl . Thus it is theoretically ex¬ 
pected that the depth of the cloud would affect the re¬ 
lation between velocity fluctuations (e.g., in the CVD) 
and projected scale. We explore this effect here using 
the simulations described in sections and 13.31 
As a limiting case of large thickness, in the statisti¬ 
cal (fractional Brownian motion) modelling mentioned 
in section [3^ we used a 128 x 128 x 128 grid (with Rel¬ 
ative Thickness TZ = z/x = 1) io generate the velocity 
field. We distribute 4000 cores randomly on the grid 
points (top panel of Figure EJ. Each core is assigned a 
centroid velocity corresponding to the velocity at its grid 
location. Thus, we assume each core perfectly traces the 


underlying velocity field. The power law index of the 
core velocities is then fitted as shown in Figure El The 
distribution of cores in the 128 x 128 x 128 grid is shown 
in the top panel of Figure El Each dot represents a core. 
The bottom panel of Figure El shows the plots of the 
power law index of the CVD-L 2 I 3 (^ 3 d) relation (7) ver¬ 
sus the power law index of the energy spectrum (/3) of 
the velocity field. Circles are results from the CVD-Z 3 D 
relation, while triangles are from the CVD-T 213 relation. 
For 1 < /3 < 3 the relation between 7 and /3 of the 3D 
scenario is close to the theoretical relation between /3 and 
7, 7 = (/3 — l)/2, however, the 2D scenario is quite dif¬ 
ferent, with the CVD-L 2 D slope being consistent with 
zero. This reflects the difference between L 2 D and 
when the thicknes s is large. For /3 > 3, 7 asymptotically 
approaches unity (|Brunt fc Heveil[2002H . 

As a limiting case of small thickness, we consider a 
128 X 128 X 4 slab (the relative thickness TZ = 1/32) 
from the previous 128 x 128 x 128 grid [TZ = 1) and put 
4000 cores randomly on the grid points (Figure fTHl) . The 
power law index of the CVD-L 2 Z 3 {ho) relation is then 
fitted based on this core sample (Figure ITOl). The distri¬ 
bution of cores in a 128 x 128 x 4 grid is shown in Fig¬ 
ure [ini The bottom panel of Figure HU] shows the power 
law index of the CVD-L 2 D{hD) relation (7) as a func¬ 
tion of the power law index of the energy spectrum {f3) 
of the velocity field. As in Figure El the circles indicate 
results from the three-dimensional CVD ((313) relation, 
while triangles are from the projected CVD relation. For 
1 < /3 < 3 the relation between 7 and (3 of both 3D and 
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Table 3 

Summary of data fitting 


Model^ 

Power law index 

Taurus cores 

0.5 ±0.05 

Taurus cores (> 6 km/s) 

0.3 ±0.02 

Taurus cores (< 6 km/s) 

0.5 ±0.3 

Taurus gas 

0.4 ±0.01 

Taurus gas 

0.5 ±0.03, 0.2 ±0.01 

Perseus cores (0 pc~ 7 pc) 

0.01 ± 0.004 

Perseus cores (7 pc-^ 13 pc) 

0.6 ±0.01 

Ophiuchus cores 

0.01 ±0.01 

RTl_L0.0625z_M14 

0.2 ±0.1 

RTl_L2.5z_M14 

0.05 ±0.11 

RTl_L5z_M14 

0.05 ±0.01 

RTl.L10z.M14 

0.1 ± 0.005 

HDl.L10z.M10 

0.05 ±0.02 


^ Models from Table 2 are only listed if a sufficient 
number of cores is identified by GAUSSCLUMPS to 
perform the analysis. 


1 8 15 23 


w 



0 5 10 15 20 

I-2D (pc) 


1 5 10 15 



Figure 5. The CVD plot of the two group of cores in Figure [4| 
top: For the case with core centroid velocity greater than 6 km/s, 
two power laws are needed to fit the data, the power law indices 
are 0.1 zb 0.03 {L 2 D < 5 pc) and 0.3 zb 0.02 (5 pc< L 2 p < 10 
pc), respectively;feo££om: the case with core centroid velocity lower 
than 6 km/s, the fitted power law index of the CVD-L 2 D relation 
is 0.5zb0.3. 

2D relation is close to the theoretical relation between P 
and 7 , 7 = (/3 — l)/2. For /3 > 3, 7 asymptotically ap¬ 
proaches unity (|Brunt &: Heved[200^ . The 3D and 2D 
scenarios are now similar. This is natural since the 2D 
length scale of a thin slab is approximately equal to the 
3D length scale. 

We studied the CVD-L 2 n relation of 
the fractional Brownian motion model 
fields of different relative thicknesses 

n = 1/32,1/16,1/10,1/8,1/5,1/4,1/2,1- The re¬ 


53 54 55 56 57 



Figure 6. The centroid velocity of cores vs. RA. in Perseus (dia¬ 
mond, with the top x axis) and Ophiuchus (circle, with the bottom 
X axis). 

suiting CyV'D-L 2 D relation with j3 = 2 are shown 
in Figure IIH with the fits to the relation and 
resulting 7 values overlayed. The panels in Fig¬ 
ure show the dependence of the /3 — 7 relation on 
the thickness of a cloud, where TZ is 1/32, 1/16, 1/12, 
1/8, 1/5, 1/4, 1/2 and 1. We find that there is a 
critical relative thickness oi TZ 1/10—1/8, above 
which the CVD-L 2 D relation significantly deviates from 
the CVD-I^d relation. This deviation is larger for 
steeper energy spectra. Since the CVD-L 2 D relation 
for Perseus and Ophiuchus is fiat below 7 pc and 3.5 
pc, their thickness/transverse extent ratio must be 
TZ > 1/10 — 1/8. With transverse extents of 16 pc 
and 8 pc for Perseus and Ophiuchus respectively, this 
corresponds to a thickness h > 1.6 pc-2.0 pc for Perseus 
and h > 0.8 pc-1.0 pc for Ophiuchus. Conversely, 
the CVD-L 2 D relation for Taurus is consistent with a 
TZ < 1/10 — 1/8 or h < 2.5 pc- 3 pc. 

The upturn of the CVD-T 2 D relation seen in 
the observations of Perseus (at 7 pc) and Ophi¬ 
uchus (at 3.5 pc) are also seen in the fractional 
Brownian motion simulation (Figure HID , albeit 
less obviously. Nonetheless, the location of the 
upturn in the fractional Brownian motion simu¬ 
lation seem to correlate with the thicknesses of 
the data cubes, which suggests that the differ- 
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Figure 7. Similar plot to Figure [T] for Perseus molecular cloud. 
The CVD data points with apparent separation 0 < L 2 D < 
7 pc can be fitted with a power law of CVD(km/s)= (1.4 ± 
+ (0.1 ± 0.04). While those points in 
7 < L 2 D ^ 13 pc can be fitted with CVD(km/s)= (0.5 ± 
0.001)1,20(pc)° ®"‘"°'°^ - (0-1 ± 0.001). The break point (7 pc) 
is determined from the lower panel by eye. 

ences in the upturn locations in the observations 
of Perseus and Ophiuchus track the differences in 
their thicknesses. 

To further constrain the depth of molecular clouds, we 
have developed a simple geometrical model of how the 
integrated linewidth of a cloud depends on its transverse 
scale and thickness. This is discussed in the next section 
(section 14.31) . 

4.3. Constraining cloud thickness from the 
linewidth-size relation 

In order to understand this critical value of the relative 
thickness TZ, we perform a simple analysis of the relation 
between the linewidth cr and 2D/3D scale L 2 D/I 3 D, and 
compare with the observational data. For a cloud with 
an energy spectrum of index /3 = 2 the integrated line- 
width is related to the 3D size by 

cr oc lli, (3) 

where I^d = h'^, h is a typical scale along the 

line of sight of the molecular cloud. So we have 



Figure 8. Similar plot to Figure[T] for Ophiuchus. The CVD data 
points with apparent separation 0 < T 213 ^ 3.5 pc can be fitted 
with a power law of CVD(km/s)= (0.8 ± 0.1)7,20 -t 

(0.05 ± 0.1). Those points in 3.5 < L 2 D ^ 5 pc can be fitted with 
CVD(km/s)= (0.5 ± 0.05)1,20(pc)° '‘*° °’' - (0.03 ± 0.001). The 
break point (3.5 pc) is determined from the lower panel by eye. 

In the limit of L 2 D 3> h, the above two relations have 
approximately the same function form. In the limit of 
L 2 D ^ h, a would not depend on L 2 d- The two func¬ 
tions are illustrated in Figure [13] The specific value of 
TZ = 1/10 — 1/8 of the threshold is not clear from this 
figure, although the general trend that the difference be¬ 
tween the a-L 2 D relation and the u-Isd relation becomes 
small at large transverse scales is clear. 

We have applied this simple geometric model to Tau¬ 
rus, Perseus, and Ophiuchus. We selected random loca¬ 
tions in Taurus and measured the average linewidth in 
groups of concentric circular regions centered on these 
locations. This analysis gives a lower limit to the overall 
thickness of Taurus. The lower envelope of the points 
on the a — L 2 d{ 13 d) plane provides a lower limit on the 
average thickness of the cloud. This lower envelope can 

be fitted with [a/vo)= {L 2 D/hY + (FigureHSj), 

and uq =0.3 km/s, h =0.7 pc. 

We applied a similar analysis to the Perseus (Figure lTd)) 
and Ophiuchus (Figure [T5|) molecular clouds and found 
that The lower envelope of the points on the (J—L 2 d( 13 d) 
plane is nearly flat in both clouds. Note that Perseus and 
Ophiuchus are elongated, so the spatial dynamic range 
for concentric circular regions average is smaller that of 


a oc 
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Figure 9. Top: Distribution of cores in a 128 x 128 x 128 grid. 
Each dot represents a core. Bottom: Relation between the power 
law index of energy spectrum (/3, see equation fT]) and the power 
law index of the CVD-L 2 D (Isd) relation ( 7 ). The 1 < /3 < 3 part 
of the solid line shows the theoretical relation between /3 and 7 , 
7 = (/3 — l)/2. For /3 > 3, 7 asymptotically approaches unity. The 
circles (3D) and triangles (2D) are the results of the model. 

CVD analysis. The fitted parameters are vq =0.6 km/s, 
h =3.6 pc for Perseus, and vq =0.4 km/s, h =2.0 pc for 
Ophiuchus. Note the thicknesses h of Perseus and Ophi- 
uchus are both lower limits, since the upturn as Taurus 
is not seen. 

In this work, the thickness of a molecular cloud has 
been estimated with the slope of the CVD-L 2 D relation, 
the fitting to the lower envelope of the a-L 2 D relation 
for all the three clouds, and the upturn location of the 
CVD-L 2 D relation for Perseus and Ophiuchus. 1) Based 
on the slope of the CVD-L 2 D relation, we conclude that 
the relative thickness TZ of Taurus is smaller than 1/10- 
1/8, while TZ of Perseus and Ophiuchus are larger than 
1/10-1/8. The relevant transverse scales are about 25 
pc, 16 pc, and 8 pc for Taurus, Perseus, and Ophiuchus, 
respectively. The estimated upper limit of the thickness 
of Taurus is 2.5 pc-3.0 pc. The estimated lower limits 
of the thickness of Perseus and Ophiuchus are 1.6 pc-2.0 
pc and 0.8 pc-1.0 pc, respectively. 2) The characteristic 
thickness h of Taurus is 0.7 pc from the fitting to the 
lower envelope of the a-L 2 D relation, which is consistent 
with the upper limit obtained in 1). The lower limits of 
the thickness h of Perseus and Ophiuchus are 3.6 pc and 
2.0 pc, respectively. 3) The thicknesses possibly indi¬ 
cated by the upturn location of the CVD-L 2 D relation of 


150 




P 


Figure 10. Similar to figure [O] for a 128 X 128 X 4 grid. 

Perseus and Ophiuchus are 7 pc and 3.5 pc, respectively. 
They are consistent with the lower limit obtained in 1) 
and 2). 

The res ults are con s istent with previous studies. For 
example, lArce et ^ (1201111 used the morphologies of 
bubbles within Perseus to conclude that its thickness is 
about 15% to 30% of its transverse scale. This is con¬ 
sistent with our conclusion that the depth of Per seus is 
large r than 15% of its projected size. In addition, (ILorenl 
1198911 find that the linewidth-size relation may not hold 
in Ophiuchus, which could partiall y explain the flat CVD 
slope we find. Enoch et al.l (1200811 find that the average 
dust extinction of Ophiuchus is much higher than that of 
other nearby clouds of similar projected size. This can 
be explained if Ophiuchus is extended along the line- 
of-sight, and thus, TZ » 1/10 — 1/8. The upturning 
point of the CVD in Figure [7] and Figure [5] may also be 
an indicator of the thickness of the clouds. In Figures 
MM and M the linewidth at the flattened part of the 
linewidth-size relation is 0.3 km/s to 0.6 km/s, which is 
larger than the thermal linewidth (at 10 K, the thermal 
linewidth of ^^CO is about 0.07 km/s). The flattening of 
the linewidth-size relation is little affected by the thermal 
broadening. 

4.4. Influence of the turbulence 

We used hydrodynamic simulations to study the im¬ 
pact of density-velocity correl ations on projection effects 
(e.g. iBrunt fc Mac Lq-^ 1200411 . The transverse scale of 
the simulated clouds is 10 pc, and we varied the cloud 
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Ljd (pc) 


Figure 11 . CVD-L 2 D plots similar to Figure [Tl for fractional 
Brownian motion model (/3 = 2, see equation^ of cores distributed 
in grids with relative thicknesses: (a)l/32; (b)l/16; (c)l/10; 
(d)l/ 8 ; (e)l/5; (f)l/4; (g)l/2; (h)l. The power law is fitted to 
the data points with 0.1 pc < L 20 < 3 pc. 

thickness from 10 pc to 0.625 pc, by considering thinner 
sub-sections of the simulated clouds. When the thickness 
is smaller than 2.5 pc {TZ < 1/4), there are not enough 
cores in the 3D data cubes, so the errors are larger. 

Figure [H] shows the CVD-L 2 D relation for cases of 
various cloud thicknesses. The CVD-L 2 _d relations with 
relative thickness TZ > 1/10 — 1/8 are nearly flat, con¬ 
sistent with the influence of thickness on the CVD-T 2 D 
relation discussed in last section. 

It is likely that the centroid velocity of cores is af- 
fected little by the global cloud Mach number (see also 
lOffner et al.ll2009all . In fact, there is a no simple corre¬ 
spondence betwee n the density distribution a nd the dis¬ 
tribution of cores (|Padoan &: Nordlundll200^ . In some 


cases, the temperature, mean density, type of driving, 
scale of driving are as important or even more impor¬ 
tant. If the CVD is primarily inherited from the cloud 
turbulent velocities, then we would not expect the CVD 
power index to depend on the global velocity dispersion, 
which simply sets the normalization of the cloud energy. 

5. CONCLUSION AND DISCUSSION 

We studied the relationship between the core velocity 
dispersion (CVD) and the projected separation of cores 
{L 2 d) in molecular clouds. We compared the results with 
an analytic model and hydrodynamic simulations. We 
come to the following conclusions. 

1) CVD is based on core centroid velocities, which 
reflects collective properties of dense gas condensa- 
tions, thus is less affect ed by density fluctuations (c.f. 
iBrunt fc Mac Lo'^ 12004(1 and more sensitive to viewing 
geometry. 

2) The relation between core velocities as a function 
of L 2 D and l^D is the same when the relative thickness 
72. of a cloud (the ratio of cloud length to cloud depth 
along the line of sight) is smaller than 1/10 to 1/8. This 
result can potentially be used to constrain the line-of- 
sight cloud dimension, which is hard to measure directly. 
A simple functional fit to the CVD index and/or h in 
equation |T] can potentially be used to provide a limit to 
the thickness (line of sight dimension) of a cloud. 

3) The CVD-L 2 d(Z 3 d) in the Taurus molecular cloud 
can be fitted by a power law with an index of 0.5, resem¬ 
bling Larson’s linewidth-size relation. This may indicate 
that the ^^CO cores in Taurus are distributed in a vol¬ 
ume that is thinner than 2.5 pc-3.0 pc, corresponding to 
1/10 to 1/8 of the transverse cloud scale. Performing the 
same analysis we conclude that the relative thickness 72 
of Perseus and Ophiuchus is larger than 1/10-1/8. The 
corresponding lower limits of the thicknesses are 1.6 pc- 
2.0 pc for Perseus, and 0.8 pc-1.0 pc for Ophiuchus. 

4) We conclude that the thickness of a molecular cloud 
can be measured by using average line profiles in concen¬ 
tric regions of different sizes, where the <j-L 2 d relation 
flattens at small scales. The characteristic thickness h of 
Taurus is 0.7 pc. The lower limits of the characteristic 
thicknesses h of Perseus and Ophiuchus are 3.6 pc and 
2.0 pc, respectively. 

Our conclusions are consistent with priors observa¬ 
tions, in which the scale and morphology of shells 
withi n the clouds ind i rectly constra ins the cloud thick¬ 
ness (jArce et al.ll2011l : iLi et al.ll2015D . This method pro¬ 
vides a promising means of probing the cloud depth in re¬ 
gions where shells arising from stellar feedback are sparse 
or absent. In future works, we will apply this analysis to 
more diverse regions, e.g. molecular clouds with different 
column densities. 
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Priority Research Program of the Chinese Academy of 
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Figure 12. Relation between the power law index of energy spectrum (^) and the power law index of the velocity field ( 7 ). The 1 < ^ < 3 
part of the solid line shows the theoretical relation between 0 and 7 , 7 = (/3 — l)/2. The circles (3D) and triangles (2D) are the results of 
the model. Relative thicknesses in each panel: (a)l/32; (b)l/16; (c)l/10; (d)l/ 8 ; (e)l/5; (f)l/4; (g)l/2; (h)l. 
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Figure 13. The solid line denotes the function (cr/i?o)= {L^ 2 D/h)‘^ + , while the dashed line represents the function ((j/vq)= 


{^ 3 D . It can be seen that the difference becomes small at large L 2 Dlh. a is normalized with a characteristic velocity vq =0.3 km/s, 
while L 2 d(^ 3 d) is normalized with thickness h =0.7 pc. The difference between these two function is plotted with dash-dotted line. The 
asterisks shows th e linewidth measured at different scales. The lower envelope of these points probes the thinnest part of the cloud (cf. 
ILi Houdell200^ . It flattens at the small scales. 



Figure 14. The same analysis of Perseus molecular cloud as Fig urell3l The characteristic velocity vq and the thickness h are 0.6 km/s 
and 3.6 pc, respectively. Since the upturn similar to that in figure fT^ is not clearly seen, here h is only a lower limit of the thickness of 
Perseus. 
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Figure 15. The same analysis of Ophiuchus molecular cloud as Fig ure 1131 The characteristic velocity uq the thickness h are 0.4 
km/s and 2.0 pc, respectively. Since the upturn similar to that in figure fT^ is not clearly seen, here h is only a lower limit of the thickness 
of Ophiuchus. 
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Figure 16. The CVD plot for a simulated cloud with a line of sight thickness of (a) 0.625 pc (with 'JZ = l/16);(b) 2.5 pc (with TZ = 1/4); 
(c) 5 pc (with 'IZ = 1/2); (d) 10 pc (with 'IZ = 1). For all panels, the Mach number is 14. The fitting range is 1 pc to 10 pc, corresponding 
to the fiat part of the energy spectrum of the hydrodynamical simulation (see figure (2]). 














